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Abstract 

The shear viscosity of a gluon gas is calculated using the Green-Kubo relation. Time correlations 
of the energy-momentum tensor in thermal equilibrium are extracted from microscopic simulations 
using a parton cascade solving various Boltzmann collision processes. We find that the pQCD 
based gluon bremsstrahlung described by Gunion-Bertsch processes significantly lowers the shear 
viscosity by a factor of 3 — 8 compared to elastic scatterings. The shear viscosity scales with the 
coupling as n ~ l/(a 2 log(l/a s )). For constant a s the shear viscosity to entropy density ratio 
n/s has no dependence on temperature. Replacing the pQCD-based collision angle distribution of 
binary scatterings by an isotropic form decreases the shear viscosity by a factor of 3. 



I. INTRODUCTION 



Results from the Relativistic Heavy Ion Collider (RHIC) at Brookhaven and the Large 
Hadron Collider (LHC) at CERN suggest the formation of a new state of matter, the quark 
gluon plasma (QGP) in ultra-relativistic heavy ion collisions. Large values of the elliptic 
flow v 2 observed in these experiments [IHZ1 lead to the indication that the QGP behaves like 
a nearly perfect fluid, i.e. its shear viscosity to entropy density ratio is small. This makes 
dissipative hydrodynamics a promising candidate to describe the collective behavior of QGP 
[8HT3]. However, hydrodynamics only works if a system is close to thermal equilibrium, 
which is not applicable for the state of matter shortly after the initial heavy-ion collision 
where transport theory based models are a more suitable approach [HHT7]. One of today's 
challenge in theoretical nuclear physics is to find a transition between the transport and 
hydrodynamic description as well as an appropriate set of transport parameters, for example 
the shear viscosity rj of the collective medium. In [18] a full leading order pQCD evaluation 
of shear viscosity in the collinear approximation of gluon splitting is reported. Calculations 
with AdS / CFT correspondence in TV" = 4 super-symmetric Yang- Mills theory in large N and 
strong coupling limit show a small universal but no n- vanishing lower bound of rj/s = 1/47T 
[TT?] . A recent analysis within a quasi-particle description and employing a Boltzmann-Vlasov 
ansatz in the relaxation time approximation also gives a small number when extrapolating 
in the non perturbative regime [20J. Considering the fully microscopic transport description 
BAMPS (Boltzmann Approach of Multi-Parton Scattering) while employing binary elastic 
99 ~~ 99 an d also inelastic bremsstrahlung gg -H- ggg pQCD collisions, the shear viscosity 
was extracted by using approaches motivated from first and second order hydrodynamics 
[2HI22]- A small viscosity to entropy ratio close to the Ads/CFT conjecture resulted, which 
is due to the incorporation of the bremsstrahlung processes being more efficient because 
of larger momentum deflection. Similar ideas have been given in [231 124] by also using 
pQCD elastic and inelastic matrix elements and modeling the equilibrium deviations of the 
distribution function in consistency with kinetic theory. Most recently, in [25] the shear 
viscosity is extracted also from a microscopic transport description employing using the 
Green-Kubo relation and perturbative QCD matrix elements for elastic scattering. 

In this work we will extract the shear viscosity of a gluon gas from microscopic transport 
calculations with BAMPS employing also directly the Green-Kubo relation. The Green- 
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Kubo relation does not rely on model assumptions or the form of equilibrium deviations. 
Thus this work can provide an independent cross check to the previously published works 
[2TJ[22]. In principle, the present calculations are rather close to those of [23]. In Sec. [lT| we 



give a brief overview of the method. In Sec. Ill the parton cascade BAMPS and numerical 



setups are briefly introduced. We demonstrate in Sec. |IV| how to extract the shear viscosity 
numerically. The results are presented in Sec. [V] and are compared to previously published 
calculations. The dependence of the shear viscosity on the pQCD coupling constant a s , 
temperature, and the distribution of the collision angle are investigated and discussed. A 
comparison with recent results presented in [25] are be discussed in Sec. V C 



II. GREEN-KUBO RELATIONS 



Green and Kubo showed in serial papers [251 [27] that transport coefficients like heat 
conductivity, shear- and bulk viscosity can be related to the correlation functions of the 
corresponding flux or tensor in thermal equilibrium. The physical motivation is given by 
Onsager's regression hypothesis [25], stating that fluctuations are present in every equili- 
brated system. Dissipation of fluctuations has the same origin as the relaxation towards 
equilibrium once the system is disturbed by an external force. Both the dissipation and the 
relaxation time scales are determined by the same transport coefficients. 

In this paper we concentrate on shear viscosity r]. The corresponding Green-Kubo relation 
has the following form [29J: 

i? = dt J v d3r <^MK J ( >°)>, (i) 

where T is the temperature, 7T 4 - 7 , i, j — x,y,z are components of the shear stress tensor, and 
(• • • ) denotes the ensemble average in thermal equilibrium. The correlation between shear 
components at time t — and at t is the sum over i and j. Furthermore one can show that 
(7r^'(r,f)7r^(0,0)) = 10(7r^(r, t)vr^(0, 0)). 

The Green-Kubo relation Q is the long wave limit k, tu — > of the Fourier transform of 
the integrated correlation function [29J . Its derivation can be found elsewhere [301432] . 

Within kinetic theory, the correlation of shear components can be computed by employing 
transport simulations for a particle system embedded in a static volume. Simulations using 
the hadron cascade model UrQMD [331 El] and the parton cascade model PCM [35] have 
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been performed to calculate the shear viscosity of a hadron gas [36] and of a quark-gluon 
plasma [25J. 

In this paper we use the parton cascade BAMPS [HI EZ] to calculate the shear viscosity 
of a gluon gas including perturbation QCD based bremsstrahlung processes gg -H- ggg. 



III. THE PARTON CASCADE BAMPS AND NUMERICAL SETUPS 

The parton cascade BAMPS solves the ultrarelativistic Boltzmann equation 

jfd^(x,p) = C(x,p) (2) 

for on-shell particles. Their interactions are simulated via Monte Carlo techniques based on 
the stochastic interpretation of transition rates [H] . Collision probabilities of two particles 
in a spatial cell of a volume of V ce \\ and within a time step At are 

D _ <?22,23 At 

^22,23 - ^relTTT 77 > {.<*) 

JVtest Vcell 

where 022,23 are the total cross section for a 2 — > 2 and a 2 — >• 3 collision, respectively, 
^rei = (pi + P2) 2 / { e 2EiE2) denotes the relative velocity of the two incoming particles with 
four momenta pi,p 2 . The probability for the back reaction 3 — > 2 is accordingly 

1 7 32 At 

32= SE^E 3 N^ t V^^ (4) 

where J 32 is a quantity corresponding to a cross section for 3 — > 2 processes and is given by 
integrating over the final states in an interaction process [TJ] ■ 

Between the collisions the particles propagate via free-streaming. Simulations of the 
space-time evolution of the particles are performed in a static box with periodic boundary 
conditions. 

In order to ensure an accurate solution of the Boltzmann equation, the volume of spatial 
cells has to be chosen as small as possible, which obviously enhances statistical fluctuations 
in the simulation. To reduce these fluctuations a test particle method [H] is introduced: the 
particle number is artificially increased by a factor of N test . Thus, the collision probability 
has to be reduced by the same factor to keep the particle mean-free path independent of the 
value of A/test. 



For the present simulations, the length of cubic cells is set to be a factor of 3 smaller 
than the mean free path of particles, and the box length is a factor of 8 larger than the cell 
length. A^ test is chosen to have on average 12 test particles in each cell. 

However, the ratio of mean free path to box length does not influence the volume averaged 
correlation function as long as particles itself stay uncorrelated, which is the case when the 
total number of particles is large enough [38J . 



IV. EXTRACTION OF SHEAR VISCOSITY 

We calculate (ir xy (r, t)n xy (0, 0)) numerically according to the Green-Kubo relation ([JJ). 
The shear component Tr xy is defined as 

7r-»(r, t) = T^(r, t) = J pV/(r, t; p) (5) 

In the numerical simulation the volume averaged shear tensor is used 

^W^ElT' ( 6 ) 

i=l 1 

where the sum is over all particles in the box at time t. V is the volume of the box. The 
correlation of Tt xy is obtained by time and ensemble average in the limit T t — > oo 

(7) 

where N t = T t /At and i = t/At. 

With (g) and the receive the following Green-Kubo relation: 

7] = / dt (7r xy (t)ir xy (0)) (8) 



10T J0 

Figure [l] shows the fluctuation of ir xy in one event (left panel) and its correlations over 
time and ensemble average (right panel) calculated using BAMPS. Results are obtained for 
elastic binary collisions assuming constant cross sections with isotropic distribution of the 
collision angle. The initial distribution is the Boltzmann distribution at a temperature of 
T = 400 MeV 

/(P) = d G e- E ' T , (9) 




time [fm/c] time [fm/c\ 



FIG. 1: (Color online) Left: Example for the equilibrium fluctuation of the 
shear component ir xy . Right: Correlations for different isotropic and 
constant cross sections. Results are normalised by ((7t xy ) 2 ). 



with the degeneracy factor da = 16 for of gluons. We neglect Bose enhancement, which is 
a quantum effect of gluons and has not been implemented in BAMPS yet. 

From Fig. [I] we see that the correlation function is an exponentially decreasing function as 
expected for solutions of the Boltzmann equation with stochastic interpreted cross-sections 
[39], 

(Tt xy {t) 7f^(0)) = (vr^(O) 7r^(0))e-* /T , (10) 

with a relaxation time scale r. The initial variance of n xy can be calculated analytically and 
is given by 

(7f"(0) 7T*"(0)> = ~ , (11) 



where e = 3g?gT 4 /tt 2 is the gluon energy density. Inserting Eqs. (10) and (11) into Eq. (JTJ 
we obtain 

V = ^er. (12) 



Equation ( 12 ) is exactly the same as derived in Refs. [2TJ HD] , if one interprets the relaxation 
time as the inverse of the transport collision rate. In general, the relaxation time scale r 
depends on temperature and individual matrix elements of interactions. Its relationship 
with microscopic scales will be shown later explicitly. The numerical task is now to find r 
by an appropriate fit to the numerically calculated correlation function. 



V. NUMERICAL RESULTS 



In this section we present results of shear viscosity. At first we assume constant cross 
sections for binary elastic collisions with isotropic distribution of the collision angle. In this 
case we are able to cross-check our numerical results with analytical calculations. Then we 
calculate the shear viscosity of gluons including pQCD-based elastic and bremsstrahlung 
processes. 



A. Isotropic cross sections 

For binary elastic collisions with isotropic distribution of the collision angle the shear 
viscosity of an ultrarelativistic Maxwell-Botzmann gas is known from the derivation in the 
Navier-Stokes approximation [4*T| |4*2"] , which is 

r] NS = 1.2654 — , (13) 



where 022 is the total cross section for binary elastic collisions. Comparing Eq. (13) to Eq. 



(12) leads to r = 1.58/(71022) = l-58A m / p , i.e., in this case the relaxation time is 1.58 times 



the mean free path. 



Figure |2| shows our results in open circles, compared with Eq. (13) represented by a 
solid line. The standard deviations are small and displayed within circle areas. We see an 
excellent agreement between our numerical results and the analytical ones, which proves the 
applicability of our numerical method within BAMPS. 

We also calculate the shear viscosity for inelastic 2 -H- 3 process with a constant cross 
section cr 2 3 and isotropic distributions of the collision angles. In this case we have [Hj 

h2 = -r-TT CT 2 3 (14) 

a G 

for computing the interaction probability P32 for the 3 — > 2 process. 

In Fig. [2] the results of the shear viscosity for such inelastic processes, 7723, are given by 
the open squares. Compared with the shear viscosity for binary elastic collisions, 7722, we 
find that r/23 is a factor of 1.5 smaller than r/22 for 022 = 023, as indicated by the dashed line. 
This finding agrees with previously reported numerical jl3j and recent analytical calculation 
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FIG. 2: (Color online) Shear viscosity for a gas with constant and isotropic 
cross sections. The open circles (squares) are results for elastic 2—7-2 
(inelastic 2 -H- 3) process. Temperature is fixed at T = 400MeV. 



B. pQCD based cross sections 

We now consider elastic and bremsstrahlung processes of gluons based on perturbation 
QCD (pQCD). The cross sections and matrix elements are the same as used in previous 
studies H3H2HI22]. 

For elastic interactions of gluons, gg -H- gg, we use the Debye screened cross section in 
small angle approximation 

d °gg^gg _ 1 

dq\ ~ 2 (ql + mlf 

The Debye screening mass is computed dynamically from the local particle distribution 

/ = f(p,x,t) via 

/d^p 1 

where dc = 16 is the gluon degeneracy factor for N c = 3. 
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(15) 



Inelastic gg -H- ggg processes are treated via an effective matrix element based on the 
work by Gunion and Bertsch [45J. Detailed balance between gluon multiplication and an- 
nihilation processes is ensured by the relation \M gg ^ggg\ 2 = d G |.M gg£M . gg | 2 . For the case of 
bremsstrahlung-like processes the matrix element employed in BAMPS reads 



2 72n 2 a 2 s 2 48vra s q^_ 

(ql + mlf ki [(k ± - q ± ) 2 + m%] 



\M 99 ^ 9 / = S , a ,,r„ 91 © (A, " r) (17) 



q_L and denote the perpendicular components of the momentum transfer and of the 
radiated gluon momentum in the center of momentum (CM) frame of the colliding particles, 
respectively 

When considering bremsstrahlung processes the LPM-effect [ISj, a coherence effect named 
after Landau, Pomeranchuk and Migdal, needs to be taken into account that leads to a 
suppression of the emission rate for low p± particles. Since such an interference effect cannot 
be incorporated directly into a semiclassical microscopic transport model such as BAMPS, 



we choose an effective approach by introducing the Theta function in (17). This implies 
that the formation time of the emitted gluon must not exceed the mean free path of the 
parent gluon g , ensuring that successive radiative processes are independent of each other. 

We calculate the shear viscosity for a temperature of T = 400 MeV. Taking the coupling 
a s as constant, the shear viscosity to the entropy density ratio rj/s will not depend on 



temperature for the considered gluon interactions in our case [21] (see later in Sec. VC). 
From next-to-leading-log calculation the shear viscosity should roughly scale with rj NNL ~ 
g~ 4 T 3 (though it has an additional logarithmic dependence on T due to Debye screening 
mass) [12] while the entropy of a quark gluon plasma should scale in leading order with s ~ 
T 3 [47] , which implies a very weak temperature dependence for T = 300MeV . . . 600MeV 

Figure [3] and Table [j] present the 77/s ratio of a gluon gas with the pQCD-based interac- 
tions, which are the main results of this paper. Because gluons are considered as Boltzmann 
particles, we have s = 4n = AdcT^/ir 2 . 

For pure gluon elastic interactions the fit function (dashed curve) agrees with the an- 
alytical result of [HI H9] at small a s . Other numerical results from [20] confirm this 
(a 2 . log(l/of a )) 1 scaling. However, in [50H52] a different scaling 77 ~ (a s log(l/a s )) _1 was 
suggested, which is a clear difference to the above findings. This different behavior stems 
from the ansatz that the collision width r _1 ~ a s Tlog(l/o; s ) is employed in their treatment 
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FIG. 3: (Color online) Shear viscosity to entropy density ratio for a gluon 
gas with pQCD-based interactions. 
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of a simplified Boltzmann collision process in the relaxation time approximation. 

If the gluon bremsstrahlung gg -H- ggg is included, the shear viscosity decreases by a 
factor of 7-8 for a s > 0.1 and a factor of 3 — 6 for a s < 0.1. The ratio of 7722 A/22+23 
is much larger than that of [18]. The reason was already pointed out in [21]. While in 
the a s — > limit the collinear elastic processes are clearly dominant, for non- vanishing a s 
the gluon radiation prefers rather large angles in the Bethe-Heitler regime [S3]. In gluon 
radiation processes the distribution of the radiation angle is close to isotropic [23 ES]- For 
the 7]/s ratio implementing both elastic and inelastic processes the fit function is found to 
be 0.0099 {a 2 s log(1.5/a s )) 1 and is shown by the solid curve in Fig. [3J This scaling holds 
for rather small values of a s . In the limit of a s — > where bremsstrahlung becomes nearly 
collinear a scaling (a^ log(l/a s )) 1 is expected for the inelastic contributions [2"T] . 

Figure [4] shows comparisons of the present results with those previously reported by 
some of the authors. The dotted curve is taken from [21] and is obtained by using a gradient 
expansion in the linearized Boltzmann equation. The deviation from thermal equilibrium 
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-(gg -B- gg) 


-(gg f-> gg Sz gg f-> 000) 


0.01 


215 ±19 


60 ±5 


0.03 


36.4 ± 1.5 


6.9 ±0.3 


0.05 


17.04 ±0.73 


2.7 ±0.1 


0.1 


5.51 ±0.17 


0.795 ± 0.025 


0.2 


2.21 ±0.07 


0.28 ±0.04 


0.3 


1.342 ±0.035 


0.166 ±0.025 


0.4 


0.952 ± 0.03 


0.114 ±0.004 


0.5 


0.74 ±0.02 


0.087 ±0.004 


0.6 


0.61 ±0.02 


0.0664 ± 0.002 



TABLE I: Data plotted in Fig. 3 



0, which is needed to calculate the shear viscosity, is approximated by ~ (x/T)(p1/E), 
where x is a constant. Actually, \ could be a function of momentum as indicated in the 
Grad's ansatz, ~ Cqit^p^p" , where Cq is a function of temperature. It seems that the 
simple choice of in [21] makes the 7]/s values lower than the present results at small a s . 
From Fig. [4] we see that the present results have a perfect agreement with those from [22] 
(dashed curve). In [22] the shear viscosity is derived by identifying the entropy production in 
hydrodynamic (Israel-Stewart) and kinetic (Boltzmann) approach. For the deviation from 
thermal equilibrium the Grad's ansatz is used in [22]. In Fig. [4] we also show new results 
(star symbols) from another work done by Reining et al. [31]. There the shear viscosity is 
calculated as a ratio of the shear pressure to the velocity gradient in a stationary shear flow 
pattern. These new results agree well with the present ones. In summary, as noted earlier 
in [21], pQCD inelastic interactions can drive the gluon matter to a strongly coupled system 
with an rj/s ratio as small as the lower bound from the AdS-CFT conjecture. 

C. Temperature and collision angle dependence of the shear viscosity 

In this subsection we study the dependence of the shear viscosity on temperature and 
the distribution of the collision angle. We only consider elastic gluon collisions in order to 
repeat the calculations performed in [25]. Our results are shown in Fig. [5] and show some 
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FIG. 4: (Color online) Comparisons of the results in Fig. [3] with those of 

Refs. EH 



discrepancies with those of [25J at lower temperatures. At temperatures above 400 MeV our 
results are rather close to theirs. 

As already done in the calculations presented in the previous subsection, the cross section 



of elastic gluon interactions is taken in the small-angle approximation (15). The 7]/s ratios 
for the elastic gluon interactions at a s = 0.3 and at three chosen temperatures are shown 
in Fig. [5] by the circles. They agree well with a constant indicated by the line. There is 
no temperature dependence of the shear viscosity to entropy ration 77/s for the considered 
interactions, which is indicated in the findings of [25]. The authors of [25] employed the same 
cross section Eq. (15) and calculated the r\j s ratio according to the Green-Kubo relation ([TJ 
by using a transport model, the BMS implementation of PCM [35]. They found an increase 
of rj/s for decreasing temperature. 



From Eq. (12) the shear viscosity r\ is proportional to the product of the energy density 
e and the relaxation time r. It is clear that the relaxation time should be a product of the 
mean free path X m f p and a function of the distribution of the collision angle F[d(9)], i.e., 
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FIG. 5: (Color online) Temperature dependence of the rj/s ratio. PQCD 
(circles) and isotropic (squares) distribution of the collision angle in elastic 
gluon interactions are taken for comparisons. 
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m/p J 



From Eq. (15) we obtain 
9ira 2 Q 



1 



a 



gg^gg 



d{6) 



2 m 2 D (l + Am 2 D /s) 
T gg^rgg 



1 do nn ^ nn 4 sin(26 l )m| ) (s + 4m 2 



D) 



a 



(18) 
(19) 



gg-^gg 



dO (ssm 2 e + Am 2 D ) 

Because momenta of gluons are distributed with f(p) ~ e~ E ^ T , it is obvious that s scales 
with T 2 , Ugg^gg with l/T 2 , and d(9) has no dependence on T. The mean free path \ m f p = 
l/n{a gg ^ tg g) is then proportional to l/T, and 77 ~ e\ m f p F[d(9)] is proportional to T 3 , like 



s. Thus, rj/s does not depend on the temperature for the case Eq. (15). 

We now turn to discuss the dependence of the shear viscosity on the distribution of the 



collision angle. In analogy to the procedure in [25] we use the total cross section (18) to 
determine the collision probability ^ and then replace the angular distribution (19) by an 
isotropic form d iso (Q) ~ sin#. The results on rj/s are also shown in Fig. |5} The 77/s ratio is 
smaller by a factor of 3 if the pQCD angular distribution is replaced by an isotropic one. 

To understand the dependence of the shear viscosity on the angular distribution, we 
calculate the transport cross section, a tr = J d6 sin 2 9da/d9, which is roughly inversely 
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proportional to rj. For the isotropic angular distribution we have a tr = (2/3)er, while for the 
pQCD angular distribution 

tr 

= a(l + a) [ln(l + l/o) - 1/(1 + o)] (20) 

°gg^gg 

with a = Am 2 D /s. On the average (s) = 18T 2 and thus (a) = 16a s /(37r) 0.5 for a s = 0.3. 
We obtain a l gg ^ gg = 0.32a gg ^ gg , which is a factor of 2 smaller than that using the isotropic 
angular distribution. By this analytical estimation the shear viscosity effectively decreases 
by a factor of two if one replaces the forward-backward peaked pQCD angular distribution 
by an isotropic one. From numerical calculations we see a factor of 3 as seen in Fig. |5j 
This situation for elastic scatterings also indicates the efficiency of gluon bremsstrahlung 
99 999 m thermalization and lowering the shear viscosity if the fact of the finite radiation 
angle at moderate a s is taken into account. 

VI. SUMMARY 

In this work we calculated the shear viscosity with the Green-Kubo relation. The corre- 
lation function is extracted from simulations by the parton cascade BAMPS. To cross-check 
the numerical method, the shear viscosity for interactions with constant cross section and 
isotropic angular distribution is calculated and the results agree perfectly with the analytical 
ones. The main results are the values olrj/s ratio of a gluon gas including a2f>3 pQCD- 
based gluon bremsstrahlung process. These results are in perfect agreement with previously 
published numbers employing other methods. In comparison to pure elastic processes the 
the rj/s ratio is lowered by a factor of 3 — 8 for a s = 0.01 — 0.6 compared to solely elastic 
collisions when gluon bremsstrahlung is included. A nice scaling of rj ~ l/[a^ ln(1.5/a s )] 
in the coupling constant a s is found as expected from standard pQCD arguments. We also 
showed that rj/s has no temperature dependence at constant a s and the change of the pQCD 
elastic interactions to those with isotropic angular distribution brings a factor of 3 decrease 
in the shear viscosity. 
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